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Abstract. One of the most important issues in disordered systems is the interplay 
of the disorder and repulsive interactions. Several recent experimental advances on 
this topic have been made with ultracold atoms, in particular the observation of 
Anderson localization, and the realization of the disordered Bose-Hubbard model. 
There are however still questions as to how to differentiate the complex insulating 
phases resulting from this interplay, and how to measure the size of the supcrfluid 
fragments that these phases entail. It has been suggested that the correlation function 
of such a system can give new insights, but so far little experimental investigation 
has been performed. Here, we show the first experimental analysis of the correlation 
function for a weakly interacting, bosonic system in a quasiperiodic lattice. We observe 
an increase in the correlation length as well as a change in shape of the correlation 
function in the derealization crossover from Anderson glass to coherent, extended 
state. In between, the experiment indicates the formation of progressively larger 
coherent fragments, consistent with a fragmented BEC, or Bose glass. 



PACS numbers: 34.50.-s, 37.10.Jk, 61.44.Fw, 67.85.Hj 
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1. Introduction 

The interplay of disorder and interactions lays at the heart of the behaviour of many 
physical systems, including superffuid helium in porous media [I], granular and thin- film 
superconductors (2HH], and light propagating in disordered media [6HE]. A central aspect 
for ultracold bosonic systems is the competition between disorder, which tends to localize 
particles, and weak repulsive interactions, which instead have a delocalizing effect. 
Wheareas disorder tends to localize non-interacting particles giving rise to Anderson 
localization j9], weak repulsive interactions can counteract this localization in order to 
minimize the energy Eventually, interactions can screen the disorder [10] and bring 
the system towards a coherent, extended ground state, i.e. a Bose-Einstein condensate 
(BEC). 

Systematic experimental studies of this interplay are difficult in condensed matter 
systems, since interactions are strong but difficult to control p], while on the other 
hand in photonic systems only non-linearities corresponding to attractive interactions 
[7J|8] have been explored in experiments. Instead, ultracold atoms in disordered 
optical potentials are a promising system for such investigations [IUH2] due to their 
unprecedented control over the disorder strength and interactions. In fact, they have 
already enabled the observation of Anderson localization for bosons in the regime of 
negligible interactions [TBI IE], and recent experiments have investigated the effect of 
interactions on the localization properties, both in the weakly interacting rj5j[16] and 
strongly correlated [TTHT9] regimes. 

Many theoretical predictions have been made about the properties of the complex 
phases appearing in these systems [20T433] . In particular, various methods to characterize 
these phases experimentally have been proposed, including measurements of transport 
properties [MJI35], condensate and/or superffuid fractions [2U [25j [27J [36] , excitation 
spectrum [36], overlap function [37], and compressibility [381 - HO] . Recent interest has 
been in the correlation properties of disordered, interacting bosonic systems [%H - H3] . 
in order to differentiate insulating phases such as the Bose glass from the superfluid 
regime. 

Here, we expand upon and extend our previous experimental work on bosons in 
a bichromatic optical lattice [16J, with an emphasis on the correlation properties of 
our system. We measure the localization properties, spatial correlations and coherence 
properties of neighbouring states as a function of the interaction energy and study the 
derealization crossover in terms of these observables. In addition, we study in detail the 
long-range decay of the correlation function of our system. Our data provide evidence of 
a change of decay behaviour at the crossover between insulating and superfluid phases, 
in agreement with theoretical predictions. 

The paper is organized as follows: In section El we give an introduction to 
the physics of a quasiperiodic lattice, and describe the expected effects of repulsive 
interactions on bosonic atoms therein. In section El we detail the experimental scheme, 
before describing the image analysis methods employed and extracted observables 
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Figure 1. Schematic of the interaction- induced derealization in a quasiperiodic 
lattice. In the non-interacting case, the eigenstates are exponentially localized for 
sufficiently large disorder, and the absolute lowest energy state is populated (a). 
Introducing very weak repulsive interactions, several of the lowest energy states are 
populated (Anderson glass, b). The energies of different states can become degenerate 
due to repulsive interactions and their shape might be modified, giving rise to the 
formation of locally coherent fragments (fragmented BEC, c), though global phase 
coherence is not restored until the entire system forms a coherent, extended state 
(BEC, d) at large interaction strengths. 



in section HJ After showing the experimental results and comparing to theoretical 
predictions in section [5J we summarize and give an outlook in section |6j 

2. Disordered phases and quasiperiodic optical lattices 

Interactions have a profound effect on disordered systems. A schematic of the effect of 
weak repulsive interactions can be seen in figure [T] for the specific case of bosons in a 
quasiperiodic lattice, as considered in this paper. Non-interacting bosons condense into 
the absolute lowest energy state of the disordered potential (figured^). The defining 
characteristic of this Anderson localized state is its exponential shape. Adding repulsive 
interactions is expected to have a delocalizing effect. This can be understood in terms of 
a screening argument [ID] . Repulsive interactions serve to smooth over the disordering 
potential in the occupied sites, providing a flatter energetic landscape on which more 
extended states can form. For very weak interactions, several low energy eigenstates of 
the non-interacting system can become populated (figure [lb)- This regime, in which 
several exponentially localized states coexist without phase coherence, is often identified 
with an Anderson glass [231 l2l]. or Lifshitz glass [26]. At larger interaction energies, 
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Figure 2. Quasi-periodic potential. The quasi-periodic potential realized in the 
experiment for lattice incommensurability /3 = 1.238... and A/ J = 6.2. The gray 
stripe shows the energy of the first band of the combined lattice, with width 2A + 4 J. 
The quasi-periodic lattice is characterized by potential wells approximately every 
— 1) ~ 4.2 lattice sites, which arise from the beating of the two lattices (grey 
dashed line). In red, the lowest energy eigenstates of the lattice potential are shown. 
Where two lattice sites are nearly equal in energy in a characteristic well, the potential 
looks locally like a double well, and there exist symmetric (red) and anti-symmetric 
(blue dashed) eigenstates, with an energy splitting of 2 J. 



an increasing number of sites is occupied, including neighbouring wells. When these 
states overlap, locally coherent fragments are expected to form (figure Hb). In this case, 
global phase coherence would not yet be restored, and the local shape of the states 
might be modified. Some authors have called this regime a 'fragmented BEC [26] or 
Bose glass [2~T j l2~? t l41] . The number of independent fragments should decrease for larger 
interaction energies, until finally, for sufficiently large interaction strengths a single, 
extended phase- coherent state is formed, that is, a macroscopic BEC (figure HH) . The 
centre of the crossover from localized to extended, coherent states is expected to occur 
when the interaction energy is comparable to the standard deviation of disorder energy. 

Our system uses a particular kind of disorder, namely a quasiperiodic potential. 
This consists of two overlapping lattices with incommensurate wavelengths. The 
resulting potential can be seen as a strong primary lattice of periodicity d = ir/ki 
which is perturbed by a weaker secondary lattice of periodicity %/k 2 {k = 27r/A, where 
A is the wavelength of the light generating the lattice). The lattice potential can then 
be written as 

V(x) = s 1 E R sin 2 (k 1 x) + s 2 (3 2 E R sm 2 (k 2 x + 0), (1) 

where E R = h 2 kf/(2M) is the recoil energy for the primary lattice (M is the atomic 
mass), (3 = k 2 /k\, and s« are the heights of the lattices in units of their recoil energies. 
The lattice spacing of such a potential is to good approximation given by that of the 
primary lattice d = Ai/2 [45]. The essential features of such a potential are visible 
in figure [2j The potential energy minima of the primary lattice are modulated by the 
second one, giving rise to characteristic wells separated on average by d/((3 — 1). The 
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additional structure given by these characteristic wells can be employed in the analysis 
of the experimental data, as described in section HI 

For non-interacting atoms, the full Hamiltonian can be mapped onto that of the 
Harper jl6] or Aubry- Andre model [4TH50] . 

H = -jJ2 (s*+i c i + c l c i+0 + A^cos(2tt^ + 0)|c/, (2) 
j j 
where j is a label for the lattice sites, and Cj give the amplitude of the Wannier state 
centered at site j. In a tight-binding model, the tunneling energy is that of the primary 
lattice, and can be calculated in terms of the experimental parameters as [51] 

J = 1.43s?- 98 exp {-2.07^} E R . (3) 

The disorder energy can be obtained from a numerical calculation as [50] 

A = 0.5s 2 /3 2 [exp (-2.18/s?' 6 )] E R . (4) 

The Aubry- Andre model displays a transition at A/J = 2 from extended to localized 
eigenstates. In an experimental realization with sufficiently large primary lattice, 
the localized regime is characterized by the absence of mobility edges, as well as 
exponentially localized eigenstates with the same localization length. It should be noted 
that this differs both from the case of a randomly disordered system, for which any non- 
zero amount of disorder is sufficient to localize the system in one dimension [52], and 
from the case of a speckle potential, for which effective mobility edges exist due to the 
correlated disorder [53J. 

The spectrum of such a quasi-periodic potential can easily be calculated and is 
shown in figure [3] for various values of the disorder strength A/J above the localization 
transition, neglecting any external confining potential. A first striking feature is the 
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appearance of minigaps in the spectrum, the lowest of which has approximately the 
same width for all values of A/ J. The appearance of minigaps can be understood from 
the distribution of lattice sites in the characteristic potential wells (see figure [2]). The 
absolute lowest energy eigenstates are those for which a lattice site coincides with the 
minimum of a potential well. In contrast, when two neighbouring lattice sites are nearly 
symmetric in potential energy in a well, the potential appears locally like a double well, 
for which the two lowest-lying eigenstates have a separation in energy of 2 J. In fact, the 
width of the lowest minigap is approximately 2 J throughout the range of A/ J shown. 
The lowest "miniband" of energies corresponds to the lowest energy eigenstates localized 
in the potential wells 4.2d apart. Since in the experiment, only the states in the first 
"miniband" are populated, we restrict our analysis to these energies and show their 
standard deviation a(E) and the extension AE of this band. The effect of a confining 
potential on the spectrum has been analysed previously in ref. [50] . 

3. Experimental methods 

In the experiment, a degenerate Bose gas of 39 K is employed in a quasiperiodic optical 
lattice. The production of a BEC of 39 K has been described in detail previously |54j. 
A broad Feshbach resonance allows a tuning of the interactions, and even a nearly 
complete cancellation [55] . In our case, a BEC of 40 000 atoms at a scattering length 
of 250ao is initially prepared in a crossed dipole trap. The condensate is loaded into 
the quasi-periodic potential while the optical trap is decompressed in about 250 ms to 
reduce the harmonic confinement, and a gravity-compensating magnetic field gradient 
is added. At the same time, the scattering length a is changed by means of the broad 
Feshbach resonance to values ranging from a < O.lao to about a = 300ao- At the end 
of this procedure, the lattice lasers give a harmonic confinement of u± = 2tt x 50 Hz 
in the radial direction. In the vertical (axial) direction, a weak confinement of 5 Hz is 
given by a weak optical trap as well as by a curvature from the gravity-compensating 
magnetic field. The primary lattice is generated by a Nd:YAG laser with a wavelength 
Ai = 1064.4 nm and a strength s± = 6.6, which is well within the tight binding regime. 
The secondary lattice is generated by a Ti:Sapphire laser of wavelength A2 = 859.6 nm 
and variable strength up to S2 = 1.2. For these experimental parameters, the separation 
of neighbouring states is given on average by d/(/3 — 1) ~ 4.2 lattice sites. 

We estimate that around 30 lattice sites, corresponding to about 7 adjacent localized 
states, are populated during the loading of the lattice. We then define a mean interaction 
energy per particle E int = gN/7 J |<^(r)| 4 d 3 r, where g = 47rh 2 a/m and (p(r) is a Gaussian 
approximation to the on-site Wannier function. We include coupling into the radial 
directions of our system, with the consequence that the interaction energy is non-linear 
in the scattering length. Though this definition of the energy is strictly valid only in the 
localized regime, comparison with a numerical simulation of our experimental procedure 
has shown that it is a good approximation for all values of the scattering length up to 
an error of 30%. 
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Figure 4. Radial temperature (a) and condensed fraction (b) as extracted from a fit 
to the radial profile. The blue circles are for A/ J = 6.2, the red squares for A/ J = 9.5. 
The grey triangles are for a single lattice (si = 5.7). The error bars denote the standard 
error of the mean. 



The loading process is adiabatic for most of the parameter range explored until 
E'int becomes sufficiently low for the system to enter the fully localized regime [56J. 
Here, several independent low-lying excited states are populated even when it would be 
energetically favourable to populate just the ground state. A loss of adiabaticity can 
be seen experimentally as a transfer of energy into the radial direction. We measure a 
radial temperature and condensed fraction by fitting the radial profiles extracted from 
the absorption images with a two-component fitting function. In our previous work [16J, 
a radial heating was seen to occur in a region of large disorder and weak interactions. In 
this work, the disorder strength is smaller and as a consequence the radial temperature 
as well as the condensate fraction are approximately constant throughout the parameter 
range explored (figure Hj) . 

4. Momentum distribution, correlation function and phase fluctuations 

The system can be characterized by analyzing its momentum distribution and derived 
Fourier transforms. These techniques are used to extract information both about the 
local shape of the wavefunction, spatial correlations, and the coherence properties of 
neighbouring states. An image of the momentum distribution is taken by absorption 
imaging with a CCD camera after 46.5 ms ballistic expansion. This time is sufficiently 
long in order to be in the "far- field" limit [57]. At the time of release, the scattering 
length is set to below Iclq in less than one ms and kept there until the Feshbach magnetic 



Correlation function of weakly interacting bosons in a disordered lattice 



8 



(a) 

o • o <0 o • o 




-3 -2 -1 1 2 3 5 10 15 20 

momentum (k ) position (lattice sites) 



Figure 5. Example of image analysis. After integration along the radial direction of 
the acquired absorption image (a), the profile of the momentum distribution is fit with 
a modulated Gaussian to recover the phase fluctuations (b). The Fourier transform 
of the square root of the profile can then be fit with three generalized exponential 
functions to extract the exponent and local length of the localized states (c). The 
correlation function g(x) is given by the FT of the momentum distribution itself, and 
can be fit with two generalized exponentials, as in (c) in order to get g(4.2d)and g(8Ad), 
or it can be fit with a generalized exponential decay up to 20 lattice sites (d). 



field is switched off 10 ms before taking the image - at this point, the system has 
expanded a sufficient amount to neglect the effect of interactions. After such a long free 
expansion without interactions, the image of the atoms that is acquired is approximately 
the in-trap momentum distribution p(k) = (^(k)^f(k)) |58j, where ^/(k) is the Fourier 
transform of the bosonic field operator &(x). In order to recover information about the 
in-trap wavefunction, we can therefore use an inverse Fourier transform. 

Due to the quasi-periodic nature of the employed lattice potential, we expect that 
for a sufficiently homogeneous system, the in-trap wavefunction can be decomposed into 
copies of the same state with real and non-negative envelope £(x), spaced by D = 4.2c?. 
The overall in-trap wavefunction can therefore be approximated as 

1>(x) = '52a j £(x-jD)e- i +', (5) 

3 

where <f>j is the local phase, and £(x) can be taken as a generalized exponential function 
exp(— \x/L\ a ). In momentum space, the magnitude of the overall wavefunction can then 
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be written as y/ p(k) = \^(k)\S(k), where 



S(k) 



E 



(6) 



is an interference term. For many envelope functions £(x), such as the generalized 
exponentials with < a < 2, the Fourier transform itself is real and non- 
negative [59], so that the inverse Fourier transform of yj p(k) can be written as £(x)aS(x). 
This is simply the convolution of the envelope of a single state £(x) with the Fourier 
transform of the interference term, S(x), which can be approximately described as 
a series of sharp peaks (approaching ^-distributions) spaced by D, with a decreasing 
amplitude and phases that depend on the local phases <pj and amplitudes aj. 

The inverse Fourier transform of the square root of the momentum distribution 
p(k) therefore allows an estimate of the average local shape of the (wave)function £(x). 
The smallest cloud size observable with the imaging system employed is about 12 pm, 
therefore there is a finite resolution also in momentum space (about fci/35). Due to 
this finite resolution, the Fourier transform has an envelope with a width of about 10 
lattice sites. This means that we can only distinguish easily up to three neighbouring 
states. The averaged wavefunction is analyzed by fitting to the sum of three generalized 
exponential functions modulated by the primary lattice 

x — jD 



/(*) 



Yl A i exp 



■ [1 + Bcosihx + 5)] 



(7) 



see figure for examples. From such a fit, the exponent a and the local extension of 
the states L s can be extracted. 

On the other hand, the inverse Fourier transform of the momentum distribution 
itself can be employed to find the correlation properties of neighbouring states. The 
momentum distribution can be related to the first order correlation function G(x', x + 
x') = (W(x')iff(x + x')). We analyze the spatially averaged correlation function, which, 
using the Wiener-Khinchin theorem, can be expressed as 

dk 



9{x) 



G(x', x + x') dx' 



2tt 



p(k)e 



(8) 



Experimentally, g(x) is recovered simply by taking the Fourier transform of the 
momentum distribution. We fit with the same generalized exponential of ([7j) and 
recover the spatially averaged correlation between states 4.2 (8.4) lattice sites apart as 
A 2 /Ai (A 3 /Ai). Also here, the finite momentum resolution limits our analysis to three 
neighbouring sites, and it follows that g(4.2d) (g(SAd)) saturates at a value around 0.6 
(0.3). 

More information about the extent and decay of the spatially averaged correlation 
function can be gained by examining the Fourier transform of the momentum 
distribution at larger distances. While the detailed structure there is not resolvable, 
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making the data there unsuited for the analysis described above, we can extract 
information about the general shape of g(x). The data is fit with a function 




exp 



x 



jD 



• [1 + Bcos(k 1 x + 5)] 



(9) 



where (3 is the correlation exponent and L g is the correlation length. This describes the 
sum of five generalized exponential functions spaced by D, with amplitudes determined 
by the shape of the correlation function. Note that this correlation function is called 
degree of coherence by Fontanesi et al. in ref. [41| . 

Finally, the effect of a fluctuating phase between neighbouring states is seen as a 
shift of the phase of the interference in the momentum distribution. We extract this 
phase by fitting the momentum distribution within the first Brillouin zone directly with 
a fitting function 



Aexp 



(k - k 



c 



2w 2 



[1 + B cos (D(k-k c ) + 



(10) 



where kc is the center of the distribution, determined by fitting the average of all images 
of a given dataset. 



5. Observed disordered regimes 

We present measurements of the quantities described in the previous section for two cuts 
of the phase diagram shown in ref. [16J. This data was taken for si = 6.6, corresponding 
to a value of the tunneling energy of J/h — 200 Hz, and for two values of the disorder 
strength, A = 6.2J and A = 9.5J. Both lattice strengths were calibrated by Bragg 
diffraction [60J, with an estimated error on A/ J of around 15%. 

The results of the analysis of localization properties as well as correlation of nearest 
and next-nearest neighbouring states are shown in figure [6j where we plot the local 
extension of the states L s , the exponent a, and the correlation function evaluated at 
4.2c? and 8Ad, as a function of the interaction energy E int . We find that for very small 
£?i n t, the states are exponentially localized, since the exponent a ~ 1, and the local 
length L s is small, consistent with the Anderson glass regime. Increasing E- m t, the local 
length increases and the exponent increases up to a > 2. Repulsive interactions therefore 
delocalize the system as expected, or alternatively, the localization crossover is shifted 
to higher values of the disorder strength A/ J when interactions are introduced into the 
system. In the localized regime, the correlation is finite but small, due to the occupation 
of independent neighbouring localized states arising from the non-adiabatic loading. As 
E int is increased, the correlation features a crossover towards larger values, signalling 
that coherence is progressively established locally over distances of first 4.2g? and then 
8.4c?. The position of this crossover is in good agreement with the prediction of the 
simple disorder screening argument, from which we expect the centre of the crossover 
to occur when E int is comparable to the standard deviation of energies in the lowest 
miniband. For A/J = 6.2 (A/ J = 9.5), this is given by 0.26 J (0.47 J). 



Correlation function of weakly interacting bosons in a disordered lattice 



11 



2.5 



3 

- 2.0 

CD 
c 
o 

g- 1.5 



1.0 



_ 1.5 
<u 

0) 

.y i.o 

TO 



0.5 



'■'"I I L 1 

-(a) . # « □ 

V * 

-U m m 

1 1 in| 1 i i 1 1 in| 1— i- 


'■'"i I 1 1 1 

c * - 

*■ ■ - 

5J CD □ ^ 

-H+rr| 1 1 1 1 llll| 1— H- 


m •s 

_• • 


(d) 

~ 

i i ■ ■ ■ 



0.6 



04 S 



0.2 

0.0 
0.6 



0.4 5 

CO- 



0.2 



0.0 



0.1 



0.1 



E ./J 



E../J 



Figure 6. Results of the Fourier analysis. The blue circles are for A/J = 6.2, the 
red squares for A/J = 9.5. The exponent a (a) and local length L (b) are extracted 
from a fit of three generalized exponential states to the FT of the square root of the 
momentum distribution. The spatially averaged correlation g(x) is extracted from the 
FT of the momentum distribution itself, and fits with three generalized exponential 
functions evalate g{x) at 4. 2d (c) and 8Ad (d). The error bars denote the standard 
error of the mean. 



Information about the phase coherence of neighbouring states can be obtained 
by measuring the phase <fi of the interference pattern in the momentum distribution 
for repeated runs of the experiment with the same parameters. If the states are 
not phase locked, changes almost randomly at each repetition of the experimental 
sequence. In figure [7] we show the standard deviation of 0, estimated from a large 
number of repetitions of the experiment. We see a decrease of the phase fluctuations 
with increasing E- mt , that nevertheless remain relatively large in the crossover region 
where the correlation increases. The fluctuations finally drop to the background value 
only when E- ia t is comparable to the full width of the lowest miniband of the non- 
interacting spectrum (dotted lines in figure [7]). These observations confirm that in the 
localized regime the states are totally independent, which together with the localization 
properties (figure [6]) indicates that the system can indeed be described as an Anderson 
glass |23,24j. The system crosses a large region of only partial coherence while becoming 
progressively less localized as E- lTLt is increased. This is consistent with the formation of 
locally coherent fragments expected for a fragmented BEC. An analogous fragmentation 
behaviour was reported in ref. [15] Ultimately, the features of a single extended, coherent 
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Figure 7. Standard deviation of the phase measured by repeating the experiment 
up to 26 times for a given set of parameters, for A/ J = 6.2 (a) and A/ J = 9.5 (b). 
The error is estimated as A<p/y/N, where N is the number of images from which the 
phase was extracted. The grey shaded bar shows the phase fluctuations measured for 
an extended system below the localization threshold. The dash-dotted line gives the 
standard deviation for a purely random distribution. The dashed (dotted) lines gives 
the standard deviation (full extension) of energies in the lowest miniband. 



state are seen, i.e. a BEC. 

Finally, we show the overall behaviour of the correlation function in figure El The 
correlation length L g increases at larger interaction energy to values larger than the 
mean separation of states (4.2c?) . It saturates at values around 6g?, consistent with the 
imaging resolution, for both quasiper iodic and single lattice potentials. The increase 
in correlation length shows that the average size of fragments found in the fragmented 
BEC regime increases with E int until presumably only a single fragment describes the 
system. 

The exponent /3 of g(x) is seen to increase from values of around 0.5 to values 
slightly larger than 1. While such an increase is qualitatively expected, the values of 
/3 are not in agreement with expectations from theory. For a three-dimensional Bose 
gas at zero temperature, we expect a transition from exponential decay (j3 = 1) in 
the insulating regime [22] to a shape of the correlation function given by the confining 
potential, > 2 in our case. In our analysis, any finite thermal component artificially 
reduces the exponent by increasing the values of the Fourier transform at small x values. 
Indeed, we observe an exponent of 1.5 or less even in the single lattice potential, for 
which the system is superfluid. In the quasiperiodic lattice, the exponent approaches 
that of the single lattice potential for large values of the interaction energy. The decrease 
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Figure 8. Decay of the correlation function, (a) Correlation length, (b) exponent 
of the generalized exponential function. The blue circles are for A/ J = 6.2, the red 
squares for A/J = 9.5. The grey triangles are for a single lattice (si = 5.7), fit to a 
single generalized exponential function. The dashed line denotes the average separation 
of states, 4. 2d The error bars denote the standard error of the mean. 



of L g and /3 at large values of E int can presumably be explained by an imperfect removal 
of interaction energy from the system during the initial stage of expansion from the 
lattice. This would lead to a broader peak in the momentum distribution, and therefore 
a narrower shape at short distances in the correlation function. 

In recent theoretical works on disordered bosonic systems, the change of shape of 
the correlation function in one dimension from exponential decay to algebraic decay was 
used to distinguish the Bose glass from the superfluid phase [^THl3] . These theoretical 
investigations have the advantage of being able to consider large system sizes, where a 
jump in the first order correlation function G(xq, x) is an indication that fragments form, 
leading to an exponential decay of g(x) in the Bose glass regime. In the experiment, 
the correlation function can only be recovered for smaller distances, due to the finite 
imaging resolution and system size. Fragments with sizes larger than approximately 2D 
cannot be distinguished from the superfluid. In this sense, the evolution of the shape 
of the correlation function can give information about the crossover from the Anderson 
glass to the fragmented BEC (where the correlation length starts to increase), but cannot 
quantify the crossover to the superfluid in our current system. Use of a higher resolution 
imaging system and eventually larger system sizes could enable the observation of the 
shape of the correlation function at larger distances. The crossover from fragmented 
BEC to superfluid could then be quantified. 
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6. Conclusions and Outlook 

In conclusion, we have characterized the entire derealization crossover of a disordered 
bosonic system caused by weak repulsive interactions through study of the spatial 
localization, phase coherence, and correlation properties. In particular, we have shown 
the first experimental determination of the shape of the correlation function in such 
a system. We find three different regimes, in agreement with theoretical predictions. 
At vanishing interaction energy, the system can be described as an Anderson glass 
(or Lifshitz glass), with exponentially localized eigenstates without phase coherence 
between them. As E- mt is increased, the local shape of the states changes, and coherence 
is gradually established, leading to an increase of the size of locally coherent fragments. 
This regime is consistent with a fragmented BEC, or Bose glass. Finally, for sufficiently 
large E int , the features of a single, extended coherent state are observed, and the system 
returns to a BEC. The position of the crossover is in good agreement with the predictions 
of a simple disorder screening argument for the lowest "miniband" . 

The techniques shown here are quite general and might be of use for further 
investigations of disordered systems. In particular, the analysis of the correlation 
function can also be used for experimental systems utilizing speckles [l3l|T5l[T9], for 
which the methods described here could easily distinguish superfluid and insulating 
phases. In the current experiment, the length scale over which the correlation function 
could be observed was primarily limited by the imaging resolution. More generally, the 
finite expansion time is expected to be a more important limitation [57], especially when 
investigating larger sized systems. The required expansion time can easily be estimated 
by considering an in-trap wavepacket with a Gaussian width of Aa;o- In momentum 
space, this corresponds to a width Ak = 1/Axq. After a ballistic expansion for time 
t exp after release of the wavepacket from the confining trap, this momentum component 
will move a distance <i ex p = ht exp /(MAxo). In order to see the features of the initial 
wavepacket, this distance must be larger than Aa;o, which implies that the expansion 
time must be larger than M(Axo) 2 /^lil- For our parameters, this is approximately 35 ms, 
less than the expansion time of 46.5 ms used. However, we must consider also our finite 
imaging resolution Ax im , which can be estimated to artificially increase the width of the 
initial cloud to (Ax ex ) 2 = (Ax ) 2 + (Aa;i m ) 2 . Given our resolution Ax im ~ 12 /im, this 
suggests a necessary expansion time of 120 ms, much longer than what is used in the 
experiment. However, a modest improvement of the imaging resolution to 5 fim would 
be sufficient to analyze the correlation function up to distances of 30 sites. 

In the future, the Fourier analysis techniques might be used to explore the regime 
of strong correlations, E- int ^> J, which can be reached by using a quasi- ID system with 
strong radial confinement. There, the Bose glass can be attributed to the cooperation 
of disorder and interactions. There is however still debate on the exact shape of 
the phase diagram in this regime, particulary concerning the possibility of reentrant 
superfluidity [T7 l[22T[24l [3T]. Furthermore, the analysis of the momentum distribution 

| This is tpp in reference [57] for a coherent wavepacket. 



Correlation function of weakly interacting bosons in a disordered lattice 



15 



using Fourier techniques would also be of use in higher dimensional systems, where it 
might be possible to use phase retrieval algorithms to reconstruct the in-trap density 
distribution in detail [6T]l62]. 
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